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Abstract 

This paper studies autonomous synchronization of k agents whose states evolve 
on SO(n), but which are only coupled through the action of their states on one 
"reference vector" in R n for each link. Thus each link conveys only partial state 
information at each time, and to reach synchronization agents must combine this 
information over time or throughout the network. A natural gradient coupling law 
for synchronization is proposed. Extensive convergence analysis of the coupled agents 
is provided, both for fixed and time- varying reference vectors. The case of SO(3) with 
fixed reference vectors is discussed in more detail. For comparison, we also treat the 
equivalent setting in IR n , i.e. with states in M n and connected agents comparing scalar 
product of their states with a reference vector. 

keywords: consensus, attitude synchronization, output coupling, special orthogonal group. 
MSC classification: 93A14, 93C10, 93B27, 93D20, 34H15, 51F25 

1 Introduction 

Synchronization and collective phenomena have recently drawn considerable attention. 
Inspired by physical and natural systems (see e.g. [29, 34, 38]), the control community 
has developed an interest in self-coordination of multi-agent systems; applications include 
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communication networks and vehicle formations (see e.g. [9, 12, 13, 14, 15, 21, 23, 25, 35]). 
The literature addresses aspects like optimal configuration, collision avoidance, nonlinear 
dynamics, interaction graph structure, distributed controller tuning, etc. 

An important distinction must be made between decentralized or distributed feedback 
control, where controllers with (sometimes implicit) reference inputs are designed for a set 
of coupled subsystems [5, 39], and autonomous coordination, where systems evolve inde- 
pendently of any external command and a collective behavior emerges from their coupling 
only. Our work belongs to the latter framework. 

A fundamental autonomous collective behavior, called "synchronization" or "consen- 
sus", is that agents connected by a restricted set of pairwise interaction links reach agree- 
ment on a common state value. This problem has been studied mainly on basis of a 
standard "linear consensus algorithm" for states in M", see e.g. [14, 26, 35]. A more gen- 
eral convexity argument is proposed by [24]. In [31], consensus is extended to nonlinear 
manifolds, motivated by problems involving e.g. data on the sphere, agents representing 
subspaces, or rigid body attitudes. Synchronization of satellite attitudes, evolving on the 
space 50(3) of rotation matrices, has attracted much attention, both in presence of a 
reference or leader [6, 17, 20, 37] and in fully cooperative (i.e. leader- and reference-less) 
framework [25, 28, 32], mostly with second-order mechanical dynamics. All these papers 
consider full state exchange 1 between interacting agents. Consensus with partial state cou- 
pling in linear systems is studied in [33, 36]: agents are coupled by linear outputs that 
are essentially equivalent to the projection of the state on the same subspace for all pairs 
of interacting agents and for all time; state agreement is reached thanks to detectable 
zero-input agent dynamics. 

The present paper introduces two innovations in this regard. First, it studies syn- 
chronization with partial state coupling on the Lie group SO(n): a basic coupling law for 
synchronization on SO(n) (see [31]) is extended to a setting where information coupling 
the agents only contains the action of their state on a "reference vector" in R n . This output 
map is inspired by [27]; however [27], to focus on noise reduction, only studies the case 
of SO (2) — that is the circle — where the state is fully determined by a scalar output 
value. Second, we consider a setting where agents' zero-input dynamics vanish, but partial 
state information involves an (a priori) different reference vector for each interaction link. 
Full state synchronization is then recovered by partial state exchanges through the whole 
network. For comparison, corresponding results for states evolving on IR n are also given. 
The main goals of the paper are 

• for practical purposes, to propose results allowing synchronization of e.g. satellite 
attitudes with reduced information transfer; 

• from a theoretical viewpoint, to illustrate implications of exchanging different pieces 
of information along different links in coordinating networks, in particular for a non- 
linear state space — the Lie group SO(n). Unlike what first intuition might suggest, 

1 0r at least full configuration exchange, for second-order systems. 
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the fact that each agent shares information about its whole state with the network 
is not sufficient for synchronization. 

Contributions are presented as follows. Section 2 formalizes the problem, proposes a 
natural cost function and derives a gradient law on SO(n) that couples the agents via 
available information. The analog problem on R n is also introduced. Section 3 starts the 
convergence analysis with 2 agents; this single-link case admits a complete analysis and 
mirrors the behavior of tracking or observer algorithms. With a varying reference vector, 
states synchronize under persistent excitation: the "averaged output map" gives access to 
full state. A single fixed reference vector is not sufficient to synchronize the states. Section 4 
presents a detailed analysis for networks of k > 2 agents. Section 4.2 gives the result for 
time- varying reference vectors, adding a technical assumption w.r.t. the 2-agent case. Local 
state synchronization for fixed reference vectors is characterized in Section 4.1. We give (i) 
a necessary condition in terms of individual agent properties, (ii) a necessary and sufficient 
condition on IR n and a sufficient condition on SO(n) in terms of the rank of a generalized 
Laplacian matrix including the effect of reference vectors, and (iii) a tighter sufficient 
condition on SO(n) which involves a matrix of size kn 2 for k agents. Several (counter- 
examples are given; in particular we show that, unlike with full state exchange, there are 
locally stable equilibria different from state synchronization even for an all-to-all interaction 
graph. Section 5 studies 50(3) with fixed reference vectors in more detail, and proposes 
an algorithm to check if a network structure satisfies a sufficient condition for output 
synchronization and state synchronization to coincide under generic reference vectors. We 
also draw attention to poor robustness of the setting where reference vectors correspond 
to relative positions of the agents in M 3 . Simulations are provided for illustration. 

1.1 Notation 

I n is the identity matrix in R nxn . X T is the matrix transpose of X G R nxm and, for 
m — n, skew(X) = |(X — X T ). 5™" 1 is the unit sphere in M. n . The Euclidean norm of 
column-vector y G M. n is denoted \\y\\2 '■= \fy T y\ The Frobenius norm of M G W ixn with 
entries m^- is ||M||jr := y^tr (M T M) = yjYlij m ij- For X £ M. nxm with column vectors 

X\,X2, • • • , x m G W l , we write rank(X) = rank(xi, x% • • • , x m ) its rank and vec(X) G M nm 
the vector obtained by stacking X\,X2, ■ ■ ■ ,x m in a single column. For a set 5 C (IR nxn ) fc 
we denote vec(5) = {vec(Xf, . . . ,X%) T : (X u . . . ,X k ) G 5} C R kn2 . Vector product of 
Xx,x% G M. 3 is Xi x X2 G M 3 . Kronecker product of matrices is denoted ®. 

SO{n) = {Q G R nxn : Q T Q = J n , det(Q) = 1} is the Lie group of n-dimensional 
rotations; so (n) = {X G M. nxn : X T = —X} is the Lie algebra of 50 (71), that is the tangent 
space to SO(n) at identity. We consider throughout the paper 50 (n) equipped with the 
standard biinvariant Riemannian metric from its embedding in IR nxn , i.e. (QQ, QQ) = 
tr(fi T 6), {QAQ6} C T Q SO(n), {Q, 6} c so(n). SO(n) acts on W l by matrix-vector 
multiplication (Q, y) y Qy. The stabilizer of y G M n with respect to this action is 
stab(y) = {Q G SO(n) : Qy = y}. For y 7^ this stabilizer is a subgroup homomorphic to 
50 (n — 1); in particular, Q G stab?/ if and only if Q T G stab?/. 
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An undirected (finite) graph Q(V,£) consists of V a finite set of vertices — in the 
following representing the agents — and £ a set of unordered vertex pairs called edges 
- representing the undirected interaction links among agents. We denote by j^S, the 
cardinality of £, and following a customary abuse of notation we write (i, j) to actually 
represent the unordered pair {i,j}- Q(V,£) can be represented by its adjacency matrix 
A = (ay) where aij = 1 if G £ and = otherwise. In particular, = for 

all i.e. A is symmetric. We consider Q without self-loops (nodes between a vertex 
and itself), so an = 0. In the following, vertices (i.e. agents) are represented by integers: 
V = {1, . . . , k}. Vertices i, j linked by edge e = G £ are denoted i = vl(e) G V and 
j = vr(e) G V, with i < j. The set of edges attached to a vertex i is denoted ed(i). The 
degree deg(z) of vertex % is the cardinality of ed(z). 

2 Problem statement and coupling law 
2.1 Partial state coupling on SO(n) 

Consider k agents with states Qi G SO(n) for i = 1,2, . . . , k, n > 2, e.g. representing 
orientations of 3-dimensional rigid bodies for n = 3. The agents evolve according to 

dlQi = Qi u i ) Z = 1, 2, . . . , , (1) 

where U{ G so(n) is agent z's control input. This left-invariant dynamics has a physical 
interpretation: control inputs are angular velocities expressed in respective body frames. 
The goal is to reach state synchronization by just appropriately coupling agents through 

Ui. 

Definition 1 State synchronization set C s is the set of states such that all agents have the 
same state, i.e. C s := {(Qi, Q 2} ■ ■ ■ , Qk) G (SO(n)) k : Q l — Q 2 — . . . — Qk}- 

The present paper considers the following restrictions on the Uj. 
(Rl) Network coupling: admissible pairwise interactions are represented by an undirected 
graph Q(V,£) such that Ui only uses information from agents j for which G £. 
(R2) Partial state coupling: interacting agents % and j exchange outputs 

Qjyij and Qjytj (2) 

with reference vectors fjij = yji G S*™ -1 given but a priori unknown. 

Rl is classical in consensus problems. Our main originality is R2; see the introduction 
for a discussion of related work. A physical interpretation of (2) is, see Fig.l: i and j 
communicate to each other the coordinates in their respective body frames of direction y^. 
One could e.g. imagine satellites comparing the position of some star on their cameras or 
how they see each other; in the latter case, y\j would be their relative position vector in 
inertial frame and Qfyij expresses it in the body frame of satellite %. 
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Definition 2 Output synchronization set C Q is the set of states such that both outputs are 
equal on each link, i.e. C Q := {(Qi,Q 2 , • • - ,Qk) £ (SO(n)) k : Qfytj = Qjyij V(i, j) G £}. 

Clearly, C s C C . The main question, investigated in this paper, is under which conditions 
output synchronization implies state synchronization. Both fixed and time-varying yij 
are considered. The challenge starts for n > 2: for n = 2, as studied in [27], 5*0(2) is 
isomorphic to the circle and (2) is equivalent to full state information. Synchronization on 
manifolds with full state coupling is analyzed e.g. in [30, 31]. 




Figure 1: Interpretation of output map (2) for SO (3). Bottom center: two agents see a star 
in direction y^-. Upper side panels: exchanged outputs = expression of yij in the respective 
body frames. Orthonormal frames (ei,e 2 ,e 3 ), (fi,f 2 ,f 3 ) indicate agent orientations. 



2.2 Gradient law on SO{n) 

A gradient law for (at least local) synchronization on SO(n) with full state coupling is 
proposed in [31]. It expresses state disagreement by cost function 

fs = \ E \\Qi-Qi\\ 2 F= E (n-te(QjQi)) (3) 

{ij)e£ (i,i)e£ 

and makes agents "rotate towards each other" according to 

|g i = -grad 0i (/ 8 ) = Q < E skew (^) ' i = l,2,...,A;. (4) 

The gradient is computed w.r.t. the product metric on SO(n) k of the biinvariant Rieman- 
nian metric for SO(n) embedded in R nxri , see Section 1.1. Algorithm (4) requires agents 
to compare relative states, QjQj. For agents only comparing outputs (2), we propose a 
gradient law which minimizes the output disagreement cost function 

fo = \j2 lltflfo-9Wl2= E (l"tr(QfM 4 ,Q,)) (5) 
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where My := yijyjj is a rank 1 projector. This yields a coupling law satisfying Rl and R2: 
|g i = -grad Qi (/ )= Q, skew(gfM ii g i ) , i = 1, 2, . . . , k . (6) 

Both dynamics (4) and (6) are right-invariant w.r.t. absolute orientation: if (Q\(t), . . . , Qk{t)) 
is a solution then (Qi(t)R, . . . , Qk{t)R) is also a solution, for any constant R G SO(n). 

Remark 1 We consider first-order dynamics for simplicity. Gradient laws are a basic 
tool that can be extended to higher-order systems, e.g. applying force/torque that is the 
gradient of a potential in second-order mechanical systems; see e.g. [30] for extending (4) 
to mechanical systems with Euler equation dynamics. 

Remark 2 Dynamics (6) can also be mapped to the space of outputs and viewed as a gradi- 
ent in the submanifold of M := ( S 171-1 ) 2#£ given by N := {(r ei , . . . , r e#£ , s ei , . . . , s e#£ ) : r e = 
Qvi(e)Ve and s e = Q T vr{e) y e for some (Qi,Q 2 , ••• ,Qk) e {SO(ri)) k , Ve G £}. The corre- 
sponding dynamics on N extend naturally to gradient dynamics of h D = ^ ee g \\r e — s e ||| 
on M for a suitable Riemannian metric. Explicitly, (6) on the output space translates to 

^r e = — skew(rj sj) r e and ^s e = — skew(sj rj) s e for e G £ . 

j£ed(vr(e)) j£ed(vl(e)) 

The appearance ofrj sj, with j ^ e, in the gradient of\\r e — s e ||| expresses inter- dependence 
of the different outputs Tj,Sj that correspond to the same state Qc this inter- dependence 
would be encoded in the geometry on N. In this interpretation, output synchronization 
corresponds to the set N f]{(r ei , . . . ,r e#e ,r ei , . . . ,r e#£ ) : r £j G S 1 ™" 1 }. Examining this 
intersection seems to be a challenging geometric problem. 

2.3 Gradient law with partial state coupling on M. n 

An analog setting with agents on M. n will be simpler to investigate and serve for comparison 
in further sections. Consider k agents with states X{ G MJ 1 for % = 1, 2, . . . , k, e.g. represent- 
ing positions of agents in W 1 , n > 1. They evolve according to J^-Xj = Ui, where Ui G W 1 
must couple the agents to reach state synchronization Xj = Xj \fi,j. The coupling is subject 
to restrictions (Rl) and (R2): the network constraint remains unchanged; partial state 
coupling is imposed by comparing, for each edge the projection of Xi and Xj on some 

vector yij G S*™" 1 C M n , i.e. communicated outputs are xjyij and xjyij. So information 
exchanged along one edge is a scalar. An extension to projections onto subspaces, with 
G M nxm , could also be considered. The natural cost function for this setting is 

9o = - ^Vij- x Jv*j) 2 = g fa ~ Xj) T My (n - Xj) 
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with Mij = yijyfj defined as in Section 2.2. The ensuing gradient law writes 

± Xi = M ij( x 3~ x i) » i = 1,2,..., A;. (7) 

One term of (7) drives the components along direction t/ij of X{ and Xj towards each other, 
leaving the rest of the state unchanged. The idea is that, by combining information along 
different directions j/y of different edges, synchronization can be achieved on the whole state 
vector. This setting essentially differs from previous work [36, 33], where state-to-output 
map is the same for all edges and state synchronization relies on zero-input dynamics 
■^Xi = Axi making an individual system detectable from its single output. 

3 Convergence for two agents 

For k = 2 agents, there is a single edge (1,2) with reference vector j/ 12 . Moreover, invariance 
w.r.t. absolute orientation (or position) allows to reduce the state of the system to the 
relative state of one agent with respect to the other one. Convergence properties on SO(n) 
and on M. n are similar: for fixed j/ 12 , synchronization appears only along the space on which 
information is projected; with varying y 12 full state synchronization is possible. 

On SO(n), defining Q = QiQ 2 leads to reduced dynamics 

full state: f t Q = 2Qskew(Q T ) (8) 
partial state: |Q = 2Q skew(Q T M l2 ) . (9) 

These expressions are also valid for a directed link, that is when e.g. Q\ follows (6) and 
Q 2 remains static. Therefore, (8) or (9) are also building blocks for tracking / observer 
algorithms: if an external pilot imposes -^Q 2 (t) = Q2 Q(t) with known input Q(t) G So(n), 
then (9) remains valid if f t Q 1 = Q 1 (skew (Q J M 12 Q 2) + fi(t) ), see also [7, 19, 22]. 

Dynamics (8) and (9) implement gradient descent for cost functions f s (Q) = \\Q T — AiIIf 
and f {Q) = \\Q T y 12 — 2/12 1 1 i respectively. In terms of Q, state and output synchronization 
map to C s = {I n } and C = {Q : f (Q) = } = stab(yi 2 ) respectively. For (8), Q — I n is 
the only stable equilibrum, see [31]. For (9) the system behaves as follows. 

Proposition 1 Consider the dynamics (9) on SO(n) with n > 2. 

(a) Let Z y = {Q : Qy i 2 = —yu} the set of maxima of f a . If y\ 2 is constant, then 
all trajectories with initial condition Q(0) G SO(n) \ Z y asymptotically converge to C Q ; 
convergence is locally exponential. Only initial points in Sq = {I n + (cos(#)-l) (2/122/12 + 
zz T ) + sin(6 ) ) (y\2Z T — zy"[ 2 ) '■ z £ S* 71-1 , z T y\ 2 = 0, 9 G (— 7r,7r)} converge to C s . 

(b) If y\ 2 (t) is a time-varying smooth function, is uniformly bounded and there exist 
T, a > such that ^ J^ +T Mi 2 (r)dr — al n is positive definite for all t > 0, then C s is the 
only asymptotically stable limit set and attracts (at least) all solutions starting in SO{n)\S\, 
where S x = {Q G SO{n) : f s (Q) > 2}. 
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Proof: (a) Take f a as Lyapunov function. Along solutions of (9), 4rf = —\\4fQ\\% = 
— (1 + y i 2 Qy 12) fo < 0. SO(ri) is compact and all functions are smooth, so by the LaSalle 
invariance principle (see e.g. [16]) all solutions converge to the set where = 0. Points 
with I+2/12Q2/12 = belong to Z y ; since they are maxima of /„, they attract no trajectories. 
So all trajectories starting outside Z y converge to f a = 0. In the neighborhood of stab (2/12), 
we have y\ 2 Qyvi > thus -j| f < — /„, which yields exponential convergence. 

A trajectory Q(t) of (9) starting at Q(0) G S'o stays in S'o for all t > 0; writing Q as 
in the definition of So, dynamics reduce to = 0, -^9 = — |sin(0). Thus 9 converges 
to for 9(0) 7^ 7r, so points in So converge to C s = {/„}. Denote So = {/„ + (cos(6*)- 
l)(yi 2 yf 2 + zz T ) + sin(0) (y 12 z T -zyf 2 ) : z G S n ~\ z T y 12 = 0, G [-7r,7r]}. Note that 
S \ S C Z r Any Q G SO(n) can be written Q = -Rt/ with i? G stab(2/i 2 ) and [/ G S ; 
this can be checked e.g. by expressing Q in an orthonormal basis with first element y± 2 and 
showing that there is at least one solution (R,U) to Q = RU for any Q G SO(n). Now 
if U*(t) is a solution of (9), then RU*(t) is also a solution of (9), for any R G stab(yi 2 ). 
Therefore solutions starting at Q(0) = RU(0) with U(0) G So converge to R; solutions 
starting at Q(0) = RU(0) with U(0) G So \ So C Z y are constant. Thus solutions starting 
outside S do not converge to C s = {/„}. 

(b) Take f s as Lyapunov function: j-J s = -2tr(f ) = 2((Qy 12 ) T (Q T y 12 ) - 1) < along 
solutions of (9), with equality at time t iff Q 2 (t) G stab (y\ 2 (t)). Points outside {I n } cannot 
be stable since I n is the only local minimum of f s , which never increases along trajectories. 

Since yi 2 (t) is time- varying, a modified version of the LaSalle invariance principle must 
be used [2]. For a given set Si, it ensures that Q(0) G SO(n) \ S± asymptotically converge 
to {/„} if, in addition to standard LaSalle conditions and technical requirements that 
automatically hold thanks to compactness of SO(n) and smoothness of (9), the following 
hold: (i) solutions starting in SO(n) \ Si at t — remain in SO(n) \ S\ for all t > 0; 
(ii) denoting by Q(t + a;Q,t) the solution at time t + a of (9) with initial condition 
Q{t) = Q, there exists for any Q G SO(n) \ (Si U {I n }) a finite time r(Q) such that 
limsup^ f s (Q(t + t(Q); Q,t),t + r(Q)) < f s (Q, t). 

Taking Si = {Q : f s (Q) > 2}, (i) holds since ^f s < along trajectories. We now show 
that (ii) holds for t(Q) = T. Assume that this is not the case, i.e. there is a Q G SO(n)\Si 
such that 

limsup f s (Q(t + a;Q,t)) = f s (Q) Vcr G [0, T] . (10) 

t— >oo 

Since ^yi2 is uniformly bounded, yi 2 (t) is uniformly continuous. Therefore -^f s (Q(t + 
a; Q,t)) is uniformly continuous (in a and t) and it is necessary for satisfying (10) that 

= limsup min 4zf 8 {Q{t + <r;Q,t)) 

t^oo <T6[0,T] 

= limsup min 2 ( y 12 (t + a) T (Q(t + cr; Q, t)) 2 y 12 (t + a) - l) . (11) 

This requires the existence, for any given e > 0, of a sequence t m — > 00 for which 1 — 
2/12 (*m + of \Q(t m + <y\Q \t m ) T ) 2 y 12 (t m + a) < e for all a G [0,T]. Note that all Q G SO(n) 
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with eigenvalues e 1 ^ , \(f)j\ < tt/2, satisfy y T Q 2 y < y T Qy for all y E S n 1 ; this is easily seen 
using the Jordan normal form of Q. Since Q ^ S±, there exists /3((Q) > such that f s (Q) = 
2(n - tr(Q)) = 2(rt - ££ =1 real(e^)) < 2 — 2/3, equivalently E" =1 cos(0 j ) > n - 1 + /3; 
this requires < ir/2 — (3 for all j. Therefore, 

£ > yi2{t m +o) T {h-Q 2 {t m +o\Qitrn))yi2{t m +<?) > ^-{yi2{t m +<y) T Q{tm+<y\Q,tm)yi2{tm.+<y)) 

for all a E [0,T] and m E N. Bounding ||Q(t + a; Q, t) - Q\\ F by f° \\iQ(l;Q,t)\ l=t+a \\ F dl 
and using standard norm properties, we then have 



< / (l-tv(M l2 (t m + a)Q(t m + a;Q,t m )))da+ / ||Q(t m + cr; Q, t m ) - Q\\ F da 
Jo Jo 

< T max(l-tr(M 12 (t m + a)Q(t m + cj;Q,t m ))) + ^ max (\\fQ(l;Q,t m )\ l=tm+(7 \\ F ) 

o-e[0,T] ct6[0,T] 

< max (TA + T 2 VA) with A = 1 - (j/i 2 (t m + cr) T Q(t m + a; Q, t m )y 12 (t m + a)) 
o-e[o,T] 

< Te + T 2 ^e. 

Thus we need = liminf^oo ^ 1 — tr(M 12 (t + a)Q)da = \immi t ^ 00 tr(B(t)(I n — Q)) 
with B(t) := ~ J T M 12 (t + cr) <icr symmetric and -B(t) — a/ n positive definite for all t > 0. 
This is only possible if Q = I n - Therefore (ii) will be satisfied. □ 

Remark 3 (i) It is still unclear whether exponential convergence to C s holds in case (b). 
(ii) The decrease of f s = f s along solutions is a particularity of k = 2. 

We briefly consider the R n case. Defining x = x 2 — X\ leads to -^x = —2x (full-state 
coupling) or -^x = —2M\ 2 x (partial-state coupling). The first system obviously converges 
to x = -vv- X\ = x 2 . Convergence with partial-state coupling is established as follows. 

• For yi2 fixed, x globally converges to (/„ — 2/122/12) x, i.e. the projection of x 2 — X\ onto 
span{j/ 12 } is driven to and the projection of X\ — x 2 onto the orthogonal complement 
of span{2/ 12 } remains constant. 

• For varying j/i 2 and under the same conditions as for SO(n), x globally converges to 
0, implying state synchronization. This standard result actually holds under weaker 
conditions, see [4]. It is currently not clear how to adapt [ !] to SO(n). This setting 
is equivalent to [36] in the special case where ^J/i 2 = SI 2/12 here, and zero-input 
individual system dynamics are J^Xj = —Qxi in [36], with fl E so(n). 
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4 Synchronization conditions for k > 2 agents 



We now turn to the general case of more than two agents. We provide an extensive analysis 
for fixed i/jj in Section 4.1 and a result for time- varying y^, similar to Proposition 1 (b), 
in Section 4.2. 

4.1 Fixed reference vectors 

We start with the case of fixed reference vectors. An immediate result of the construction 
of our gradient law is the stability of the output synchronization set. 

Proposition 2 The output synchronization set is asymptotically stable under (6). 

Proof: Since (6) is a gradient descent system for f and f is an analytic function, the set 
C of global minima of f a is asymptotically stable [1]. □ 
The same clearly holds for (7), the corresponding system in M n . However on SO(n), as 
will be illustrated in Section 5, convergence is not necessarily exponential for k > 2. 

4.1.1 Necessary condition for state synchronization 

The following theorem formalizes the natural condition: to reach state synchronization, 
each agent's state must be observable by combining information on all edges incident to 
the agent. Later we show that this is not sufficient (at all) for C s = C Q . The necessary 
condition for C s = C Q in Theorem 1 can be extended in straightforward way to the following 
necessary condition: considering any partition of V into V\ and V2, the information on edges 
connecting V\ to V2 must make a full state G SO(n) observable between the two subsets, 
i.e. there must be at least n — 1 links on the cut with linear independent reference vectors. 

Theorem 1 The output map from state space (SO(n)) k (resp. (M n ) fc j to output space 
(iS™" 1 ) 2 ^ (resp. M?# e ) is injective if and only if each i G V has at least n — 1 (resp. n) 
linearly independent yij with = 1. If injectivity is not satisfied, then for any (Qi, . . . , Qk) 
in C Q there is a continuum of states in C Q \ C s giving the same output as (Qi, . . . , Qk)- 

Proof: We first discuss the SO(n) case. Consider an agent % and number the agents 
such that Oij = 1 if and only if j < m. Denote by h t : SO{n) x ( 1 S n ~ 1 ) m ->■ (S^ 1 )™ 
the action of Q t on y t := (y a , . . . , y im ), i.e. h^Q^yi) = (Qjyn, ■ ■ ■ ,Qfy im ). The map 
Qi H- hi(Qi,yi) is the output map of agent i. Define stab/ 1 (y i ) := {Q G SO(ri) : hi(Q,yi) = 
§i} — CXjLi stab(?/y). Thus by definition yn,..., yim are eigenvectors with eigenvalue 1 of 
all Q G stabfe(|/i) C SO(n). Then stabh(yj) = {/„} if and only if {yn, . . . ,yi m } contains 
(n — 1) linearly independent vectors. Indeed, the condition implies that all Q G stabh{yi) 
have (n — 1) eigenvalues 1, and Q G SO(n) implies that the remaining eigenvalue also equals 
1; conversely, if there is a > 2-dimensional subspace Sq orthogonal to {yn, . . . , yi m }, then 
all Q representing rotations in So belong to stab/j^j). By invariance of linear independence 
under a common rotation, the same holds if i/i is replaced by hi(Q*,yi) for any Q* G SO(n). 
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Therefore a particular output hi(Q*,yi) of agent % corresponds to a unique state Q* of agent 
% if and only if {yu, . . . ,yi m } contains {n — 1) linearly independent vectors. As the whole 
output map is just the collection of hi for all i G V, with hi depending on the state of agent 
i only, this proves the first claim. As said above, if the condition does not hold for agent 
i then stabh(yi) contains rotations in a subspace So, that form a Lie subgroup of SO(n) 
of dimension > 1. If (Qi, . . . , Qk) G C , then (Q\, . . . , QQi, ■ ■ ■ , Qk) £ C Q yields the same 
output for all Q G stab(^j), but at most one Q corresponds to state synchronization. 

The analog property for M. n is straightforward: projections on n linearly independent 
reference vectors specify the state completely, while less than n projections leave the pos- 
sibility to move the state in one direction without affecting the output. □ 

The following Examples 1 and 2 illustrate that injectivity of the output map does 
not guarantee C s = C Q . This happens because there exist output values Qjyij for states 
in C a , that cannot be obtained with states in C s . Therefore Theorem 1 only gives a 
necessary condition: unlike what first intuition might suggest, each agent injecting full 
state information into the network does not imply state synchronization. This suggests 
that a theory of observability in networks is non-trivial. 

Example 1 Consider 3 agents on SO (3) with a fully connected interaction graph and 
Ui2 = (1 5 0,0) T , ?/23 = (0,1, 0) T , ?/i3 = (0,0, 1) T . Assumptions of Theorem 1 hold, so 
each output synchronization value corresponds to a unique point in state space. Let output 
synchronization be reached at Q\yyi = Q\yvi — (1,0, 0) T , Qjyn = Qjyn = (0,0, 1) T and 
Q\y2Z = — (0, —1, 0) T . The corresponding unique state is Q\ = diag(l, 1, 1) , Q 2 = 

diag(l, —1, —1) , Q 3 = diag(— 1, —1, 1) which is not state synchronization. 

Example 2 Consider 4 agents with states Xi G M 2 where S = { (1,2), (1,4), (2,3), (3,4) } 
and with y± 2 = (1,0) T , yu = (0, 1) T , 2/23 = (cos^,sin^) T ; y 34: = (cos0, smcf)) T , assuming 
9,<f) G" {kn/2 : k G Z} and y 23 7^ ±^34. The setting is illustrated in Fig. 2. Assumptions 
of Theorem 1 hold, so the output map is injective. Note that the generalized necessary 
condition for C s = C discussed before Theorem 1 is also satisfied. Nevertheless, states 
of the form %i = (a,b) T , x 2 = (a, d + f^) T ', x 3 = (c,d) T , x 4 = (c + (d — 6)tan</>, b) T 
correspond to output synchronization for any (a, b, c, d) G M 4 , while state synchronization 
only holds if c = a and d = b. 



Figure 2: Example 2 - each agent sharing its whole state with the network is not sufficient 
to ensure that output synchronization implies state synchronization. 




(a,b) T = xi 



52 = (a,d+^) T 



(c + (d — b) tan (j>, b) T = X4 



x 3 = (c, d) T 
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4.1.2 Sufficient conditions for state synchronization 

The output cost functions can be written as a generalization of the Laplacian-based 
quadratic form often used in the context of synchronization algorithms: 



f (Qi,...,Qk) 



g (xi, ■ ■ ■ ,Xk) 



( 


/Qi\ 


Y 


vec 




( 


V 


\Qk) 


/ 




T 


( Xi \ 




IP 








\Xk) 





f 


(Qi\ 


\ 


(In ® L 9 ) 


vec 










\QkJ 


J 



(12) 



(13) 



where L 9 = (Lf) 



G M. knxkn is an extended Laplacian matrix with L?- G M nxn defined by 

for i ^ j , L 9 U = X^jLi a ijMij ■ The standard Laplacian matrix L 
has the same definition but with Mij replaced by 1, i.e. Ly G M. For undirected graphs, 
L = BB T where B G W kx # £ is the incidence matrix of Q: each column corresponds 
to an edge e and contains all zeros except 1 on row vr(e) and —1 on row vl(e). It is 
straightforward to check that L 9 = (B ® I n ) W (B <g> I n ) T where W G R n (# £ ) Xn (# £ ) is block 
diagonal with Mjj in the block corresponding to edge (i, j). Standard rank properties yield 



rank(L 3 ) < min(rankW / , n rank B) = min(^£, n rankL) . 



(14) 



The rank of L 9 is thus limited by number of edges and graph connectivity: for a graph on 
k vertices with m connected components, rankL = k — m. However (14) is far from being 
tight, as it contains no information about the output map while L 9 does: if e.g. all are 
equal, then rankL 9 = rankL < [k — 1). For agents in W 1 , L 9 taking its maximal possible 
rank is necessary and sufficient for global, exponential state synchronization. 

Proposition 3 C a = C s if and only i/rank L 9 = n(k — 1) for partial state coupling on W 1 . 
Moreover, if rank L 9 = n(k — 1) the set C Q = C s is globally exponentially stable under (7). 



Proof: Rewriting (7) using (13) we get -^x = —L 9 x with x = (xj, . . . ,x^) G K . As 
L 9 is positive semidefinite, x converges to ker(L 9 ) = C Q ; the latter is stable by Proposition 
2. Since C s C C Q = ker(L 9 ) and C s has dimension n, we have C Q = C s if and only if 
rankL 9 = n(k — 1). Linearity implies that convergence is exponential and global. □ 
The situation on SO(n) is quite different: (i) the rank condition is only sufficient; (ii) 
state synchronization is usually not global, even for a complete graph satisfying the rank 
condition; and (iii) convergence is not necessarily exponential. 

Theorem 2 Consider dynamics (6) for partial state coupling on SO(ri). 

(a) If rankL 9 = n(k — 1) then C Q = C s and the set C Q = C s is locally exponentially stable. 

(b) Denote V G W 1 kx fcn ( n_1 )/ 2 a matrix whose column vectors are an orthonormal basis for 
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vec(so(n) fc ) C ]R n k . C s is locally exponentially stable if and only if rank(V T (/ n ®L 9 )V) = 

NB1: N on- exponential asymptotic stability of C s is still possible if the rank condition fails. 
NB2: We refer to Section 1.1 for the definition of vec(so(n) k ). 

Proof: Extend f (X 1 , . . . , X k ) to Xi G IR nxn . Denote S the subspace of global minima of f 
on (M nxn ) fc . Obviously C Q C S. In addition, S contains the subspace C r s := {(X, ...,X) : 
X G R nxn }. If rank(La) = n (k - 1) then rank(J n <g> L 9 ) = n 2 k - n 2 and from (12), S 
is n 2 -dimensional, like CJ, thus S 1 = CJ. Then C Q d S = C r s and the intersection with 
(SO(n)) k yields C = C s . 

We apply Lemma 1 from the appendix. For (Q, . . . , Q) G C s , we have L" = (Ik <S> 
Q T )L 9 (I/ : (E)Q), with L 55 as defined in Lemma 1. Thanks to right-invariance w.r.t. absolute 
orientation, we can reduce the investigation to Q = I n without loss of generality, so L" = 
L 9 . From the previous paragraph, if rankL 9 = n(k — 1) then I n <8> L 9 is positive definite 
on the subspace of M. kn orthogonal to vec(CJ), thus in particular on the subspace of 
vec(so(n) fc ) orthogonal to vec(CJ). So the condition of Lemma 1 is satisfied and we get 
exponential stability of C Q , proving (a). By definition, Ly := V T (I n ®L 9 )V is the restriction 
of (7 n <g> L 9 ) to vec(so(n) k ). Obviously, vec(CJ) n vec(so(n) fc ) = {vec((fi T , . . . , n T ) T ) : Q G 
so(n)} C ker(Ly). Since the dimension of vec (C r s ) n vec(so(n) fc ) is it follows that 

Ly is full rank in the subspace of vec(so(n) fc ) orthogonal to vec(CJ) if and only if Ly has 
rank kn( ^^ _ n ( n ~ 1 ) _ ^ _ 1) Tt ( n ~ 1 ) ; proving (b). □ 

Conditions (a) and (b) of Theorem 2 are complementary. In condition (a), roles of 
graph connectivity and reference vectors can be visualized to some extent. In particular, 
from (14), rankL 5 = n(k — 1) requires at least n(k — 1) edges and thus k > 2n agents. 
Moreover, checking (a) involves L 9 G M nkxnk and is less costly than (b) which involves 
n 2 k x n 2 k matrices. On the other hand, (b) examines the rank of (I n £g> L 9 ) on the tangent 
space of (SO(n)) k and therefore is tighter than (a) which basically checks it in R n kxn k . 
Condition (b) can even prove local exponential stability of C s in situations where C D ^ C s . 
If condition (b) fails, C s can still be asymptotically stable, although convergence cannot 
be exponential. Section 5 contains examples illustrating these properties for 5*0(3). 

Example 3 On SO (2) = S 1 , an output Qjyij is equivalent to full state knowledge: ex- 
plicitly writing (6) and (4) in terms of rotation angles yields the same result up to gain 
factor 2. It is in fact equivalent to Kuramoto-like dynamics = X] a «i si n (% — 9i) 
for which local (sometimes almost- global) state synchronization holds as soon as Q is con- 
nected, see e.g. [31]. There are many synchronizing situations where, for Theorem 2(a), 
rank(L 9 ) < 2(k — 1), e.g. when Q has less than 2(k — 1) edges. In contrast, Theorem 2(b) 
perfectly captures the situation on S0(2): calculations with V constructed e.g. from basis 
elements (• • • , (o o) > (-1 o) » (o o) > ■ ■ ■) °f so ( 2 ) k show that actually V T (I n ®L 9 )V = L. 
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4.1.3 About almost-global convergence 

On IR™, stability of state synchronization is always global, see Proposition 3. On SO(n) 
with full state exchange, see [31], C s is not always (almost-)globally asymptotically stable, 
but (4) has no other stable equilibria than C s when Q is e.g. all-to-all or a tree. For partial 
state coupling on SO(n), n > 2, (6) does not retain these global properties 2 for trees and 
complete graphs. For a tree, the number of links is insufficient to share all state information: 
the necessary condition from Theorem 1 is not satisfied. The following example illustrates 
that even with all-to-all interaction and conditions of Theorem 2 holding, stable equilibria 
besides the state consensus set can exist. 

Example 4 Consider k = 7 agents on SO(3), interconnected with: y^ = (0, 1, 0) T if 
{i — j) mod 7 = ±1, i/ij = (0, 0, 1) T if {i — j) mod 7 = ±3 and y^ = (sin(a), 0, cos(a)) T if 
{i — j) mod 7 = ±2, with a = 0.04. This implies rank(L 9 ) = 18 = n[k — 1), so Theorem 2 
applies and C s = C Q is locally stable. Nevertheless, (6) admits an equilibrium ^ C s of the 
form Qi = Q y Q z {i^f), i — l,...,7, with 

I cos(#) sin(0) \ / cos(/3) sin(/3) \ 

Q z (6)=[ ~sin(#) cos(#) , Q y = 10 

\ 1 / \ -sin(/3) cos(/3) / 

and j3 = —0.0558. Checking the Hessian at that point shows that it is stable. The small 
but nonzero values of a and consequently [5, are necessary to have at the same time L 9 of 
maximal rank and stability of the equilibrium. 

4.2 Time-varying reference vectors 

Time-varying y^ are favorable for state synchronization, as they give access to different 
partial state information over time. We assume smooth functions of time IR — > S*™" 1 . 
Unlike for k = 2 agents, f s does not decrease along solutions of (6) in general. This 
somewhat complicates the analysis and forces us to introduce a scaling parameter e in (6), 
i.e. 

iQi = eQi Yl skew(gfM ij (t)g i ) , i = 1, 2, . . . , k . (15) 
{r-(ij)es} 

To show convergence to C s , we need a persistent excitation condition on the outputs. 

Assumption 1 M^- := lim^oo i f M; L j{s) ds exists and is positive definite V(i, j) E S . 
Moreover 3 , there exists Z : [0, +oo) — > [0, +oo) with lim^oo \Z{t) = and such that for 
all t and V(z, j) G £ , it holds || fi +J '(M y (s) - My) ds\\ F < Z(T) for allT>0. 

2 Contrary to an erroneous statement in [18], Theorem 2(b). 

3 This condition of uniform convergence to Mij was not acknowledged in [18]. 
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Assumption 1 ensures a well-defined "average system" (for e = 1) 

iQi = Qi E skew^fM^.) (16) 

which converges like full-state-coupling case (4), see [31]. For sufficiently small e, an 
averaging argument shows that this convergence extends to time- varying dynamics (15). 

Proposition 4 If Q is connected and Assumption 1 holds, then C s is locally exponentially 
stable for (16). 

Proof: (16) is a gradient system for the average cost f(Qi, . . . , Qk) — \ ~Y%,j)es \\^ij{Qi~ 
Qj)\\%- A straightforward extension of the arguments in [31] yields that C s is locally asymp- 
totically stable, because C s is an isolated minimum set of cost function /. Exponential 
stability holds similarly to other results in the paper. □ 

Theorem 3 If Q is connected and Assumption 1 holds, then for sufficiently small s > 0, 
C s is locally exponentially stable for (15). 

Proof: We first project the system onto one where C s is represented by a single point. 
For this we use reductive homogeneous spaces, see e.g. [8]. 

C s is a closed subgroup of SO(n) h . Define the compact homogeneous space Ai = 
SO(n) h jC s , with C s acting on SO(n) k by right multiplication, and canonical projection 
7r: SO(n) k — > Ai. The image of C s under tt is a point p G Ai. Since SO(n) k is compact, 
Ai is a reductive homogeneous space. Equip Ai with the normal Riemannian metric 
induced by the product metric on SO{n) k . Denote Fso(Qi, ■ ■ ■ , Qk, t) and F so (Qi, ■ ■ ■ , Qk) 
respectively the vector fields corresponding to (15) with e — 1 and to (16). Thanks to 
invariance of the dynamics under right multiplication of SO(n) k by C s , they respectively 
induce vector fields F(x,t) and F(x) on Ai, where F is the time average of F. Moreover, 
for both time-varying and time-invariant dynamics, point p G Ai has the same stability 
properties as set C s C SO(n) k . So p is exponentially stable under F and we must prove 
that this remains true under eF(x,t) for e small. We therefore apply Theorem 3 in [3]. 

Since all necessary regularity conditions for Theorem 3 of [3] hold, it remains to examine 
|| J t * 2 F(x, t) — F(x) dt\\ for any fixed x G Ai in a neighborhood of p. 

F(x,t)-F(x)dt < / F so {Qi,...,Qk,t)-Fso{Qi,---,Qk)dt 
i Jti 

by projection from Tq 1 

■1-2 



Qk SO{n) k to T X M, with 7r(Qi, . . 
Fso{Qi, ■ ■ ■ , Qk, t) — F S o(Qi, ■ ■ ■ , Qk) 



Qk) 

2 

dt 



x. Now 



k 



E 



fc rt2 

<hj skew (Qf / (My (t) - I tJ ) dt Q j ] 
Ju 



< 



i=i j=i 



skew (Qj / (M tj (t) - Mij) dtQj 



(17) 
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Using || skew(Qf PQj) \\ F = ±||P(I n - 
||-P||f \\Qi — Qj\\f for any symmetric P G 



(In-QjQjyP\\F< \\P\\ F \\In 

:n , and Assumption 1, we get 



QjQl 



[Eq.(17)] < k 



k k 

i=i i=i 



*2 



(M«(t) 



My) (it 



WQi-Q. 



2 

3\\F 



=1 J=l 



i=i j=i 



Since HQj — QjH^ is a quadratic form on 



pnxnNfc 



; with kernel = {(X, . . .,X) : 



X G M nxri }, there is a constant c > such that, with dist( 



the Euclidean dis- 



tance, Ei, 7 =i IIQ< 



QjHf < cdistf 



. . . , Qk), C r s ). As C s C CJ and we use 



the Riemannian metric on SO(n) induced by the Euclidean metric on M nxn , we have 
dist (R nxn)fc((Qi, . . . , Qk), C r s ) < dist 50(n) fc((<3i, • • • , Qk), C s ) with dist 50(n) fe the Riemannian 



QjWf 
t 2 



< c dist 



SO(n) 



Qk), C s 



distance on SO(n) k . This yields Y^j=i WQi 
c dist 2 vt (x,p) with p = tt(C s ). Then overall, 

dist^ (x,p) . Using local charts around p, this shows that the last condition of Theorem 3 
in [3] holds, ensuring local exponential stability of p under F(x, -t) for e sufficiently small. 



£F{z,t)-F{z)dt 



< kyfcZ{t 2 - t ± ) 



We conclude by change of timescale: 



di 



x 



e F(x,t) is equivalent to 



dr 



X 



F(x, ±r) 



□ 



It is possible to give quantitative estimates for £ based on results in [3], but these 
technical aspects are not the focus of the present paper. 

Remark 4 The conditions of Theorem 3 can be relaxed to require strict positive definite- 
ness of Mij, in Assumption 1, only on the edges of a spanning tree Q s of the interaction 
graph Q. Indeed, taking only edges of Q tree , the result of Proposition 4 directly applies; in 
particular, C s is an isolated set of minima for the cost f tree of the average gradient system 
associated to Qtree and is locally exponentially stable for that average system. Adding the 
contribution of edges in Q\Q tTee , for which averaged coupling is only partial, C s remains an 
isolated minimum set for the cost f of the average system associated to Q; then it is still lo- 
cally exponentially stable under the gradient system, so that the conclusion of Proposition 4 
still holds. The proof of Theorem 3 can be repeated as such. 

It should also be possible to extend the result to networks where a majority of fixed yij 
are combined with a few time-varying that enhance the connection on particular links: 
it suffices to show that C s is locally exponentially stable for the average system, and then 
the conclusion of Theorem 3 holds. 



5 State synchronization on 50(3) 

We pursue the analysis for fixed on special case S0(3), which has practical importance 
as describing rotations in three-dimensional space. It is also the simplest non-trivial output 
consensus setting. Qjyij = Qjyij fixes QiQj up to a rotation with axis yy, so the coupling 
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on one edge contains 2-dimensional information about 3-dimensional relative errors QiQj G 
50(3). We give further conditions to have C Q = C s , study error sensitivity highlighting 
poor robustness when yij are relative positions in M n , and illustrate with simulations. 



5.1 Weaker conditions for C Q = C s 

The conditions of Theorem 2 give too little insight into respective roles of graph and 
reference vectors to be efficiently used for network design. The following conditions, weaker 
than Theorem 2(a), ensure C Q = C s on 50(3) for small k. These small networks can be 
used as building blocks for larger networks in which C Q = C s with generic reference vectors. 

Theorem 4 Consider k agents applying (6) on 50(3) with fixed y^. 

(a) If k = 3, £ = {(1,2), (1,3), (2,3)} and rank(j/i 2 , yis, y 23 ) = rank(j/i 2 , j/i 3 ) = 
r&nk(y 12 ,y 23 ) = rank(j/i 3 , y 23 ) = 2, then C Q = C s . 

(b) Ifk = 3,£ = {(1,2), (1,3), (2,3)} and rank(j/ 12 , y 13 , y 23 ) = 3, then C Q ^ C s . However, 
output synchronization still implies state synchronization locally. 

(c) Ifk = 4,SD {(1,2), (1,3), (1,4), (2,3), (2,4)}, ran% , y b , y c ) = 3 for each triple of 



different a, 6, c G £, and (m^Tfai^) , (^WW^m) then C = C S . 

JJ ' ' ' 2% (2/12* 2/13) 7 2/24^12X2/14) ' 

Proof: Define relative states Rj = QjQj, j = 2, . . . , k. State synchronization means R 
I 3 , j = 2, . . . , k. For (a) and (b), output synchronization means R 2 yi 2 = 2/12, -R3J/13 = 2/13 
and R 3 R 2 y 23 = 2/23- For (c), output synchronization requires these on subnetwork A of 
agents {1, 2, 3} and the same conditions on subnetwork B of agents {1, 2, 4}. 

(a) Output synchronization implies yf 2 R 3 y 23 = 2/12-^3^3-^2/23 = 2/122/23- The rank 
conditions yield y 23 = ay u +byi 3 , with a,b e R\{0}. Thus yl 2 Rjy 23 -2/122/23 = a(yJ 2 R 3 yi 2 - 
1) = which implies R 3 y\ 2 = y\ 2 . Then R 3 has two linearly independent eigenvectors for 
the eigenvalue 1, hence R 3 = I 3 (see proof of Theorem 1). A similar argument yields 
R2 = h- 

(b) For any x G 5 2 we have {Rx : R G stab(yi 2 )} = {y G 5 2 : y\ 2 y = y\ 2 x}, and for 
x ^ ±2/12 there is a bijection between R and y. The rank conditions allow to write y 23 = 
aiei + a 2 e 2 + a 3 e 3 in orthonormal frame (ei,e 2 ,eq) := (?/i 2 , yi2Xyi: \ yi ^ x (v^xyi3) \ w jth 

00 \ 11 i! vwi^; ||yi2Xyi3||a ' \\y12xy13h n 

a 2 7^ 0. Define y* = — a 2 e 2 + 0363. Since y^y* = 2/122/23 7^ ilj there is a unique R 2 * G 
stab(2/i 2 ) \ ^3 such that R^yiz = 2/*- Similarly, using yj 3 {y\ 2 x y 13 ) = 0, there is a unique 
R 3 * G stab(j/i 3 ) \ I 3 such that -R 3 *2/ 23 = y*. Then R 3 *R^y 23 = y 23l so i? 2 * and #3* differ 
from ^3 but satisfy all conditions of output synchronization, i.e. C Q \C S ^ 0. Furthermore, 
states in C Q \ C s all correspond to the same relative state (R 2 *,R 3 *). Indeed, assume 
y Q = R 2o y 23 = Rj Q y 23 for some (R 2o ,R 3o ) G stab(j/i 2 ) x stabQ/i 3 ), characterizing output 
synchronization. The linearly independent requirements y\ 2 y Q = y\ 2 y 23 and y\ 3 y Q = 2/132/23 
have two solutions in 5 2 : y 23 and y*. As each admissible y Q maps to unique (R 2o ,R 3o ) G 
stab(yi 2 ) x stab(yi 3 ), the only possible (R 2o ,R 3o ) are (I 3 ,I 3 ) and (R 2 *,R 3 *). Thus C Q 
consists of two separated sets: C s and C D \C S = {{Q, R 2 *Q, R 3 *Q) '■ Q G 50(3)}. 

(c) From case (b), output synchronization in subnetwork A = {1, 2, 3} requires (R 2 , R 3 ) G 
{(R 2 *, R 3 *), (I 3 , 1 3 )}, with I 3 {R 2 *, R 3 *}. A similar argument yields constraint (R 2 , i? 4 ) G 
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{(i?2*, Ra*), (I3, h)}, with ^3 ^ {i?2*,-R4*}, for output synchronization in subnetwork 
B = {1,2,4}; hereafter we characterize {R 2 +,Ra*) by y+ G S 2 . Output synchronization 
in the whole network combines these conditions. Then consistency for R 2 requires either 
R2 = h (implying R 3 = R4 = I 3 , state synchronization) or R2 = R 2 + = R2*- We show 
that the latter is impossible. 

Let (di, a 2 , 03) coordinates of 2/23 as in (b), in particular a 2 7^ 0. Write y 2i = t>i e i+b 2 e^+b 3 e 5 
in orthonormal frame (ei,e 4 ,e0 = (yi2, u yi2XVl4 u , yi jf x ^ 2Xy 14 > ) with b 2 ^ by the same 
argument as in (b). Define i? G stab(yi2) with (ei,e4, es) = i? (ei, e2, 63). Note that 
RL(V23 ~ oiei) = (y, - aiei) = (a 2 e 2 + a 3 e 3 ) and R T Rl+{y 2i ~ Mi) = R T {y+ ~ Mi) = 
(6 2 e 2 + b 3 e 3 ). Assume R 2 + = R 2 *. Then R^(y 23 ~ ai^i) x R T R^(y 2i - Mi) 

= (a 2 e 2 + a 3 e 3 ) x (6 2 e 2 + b 3 e 3 ) = mei = me x = ( (a 2 e 2 + a 3 e 3 ) x (6 2 e 2 + b 3 e 3 ) ) 

= (i^(a 2 e 2 + a 3 e 3 )) x (R^(b 2 e 2 + 6 3 e 3 )) = (-a 2 e 2 + a 3 e 3 ) x (-b 2 e 2 + b 3 e 3 ) . 

The first and last expressions yield ^ = g. Expressing 02, 03, 62; ^3 by scalar products leads 
to fai2xy 23 ) r fai 2 xy 13 ) = (y 12 xjQ4) T (yi2xyi4) violating the last assump tion of the statement. 

Thus i?2* = -^2* is impossible and C D \ C s = 0. □ 

Remark 5 (%) i?an& conditions in Theorem 4(b) and (c) are only violated on a set of mea- 
sure zero for freely chosen yij G S 2 . Situation (a) — coplanar but non-aligned yi 2 ,yi 3 ,y 23 
- is only violated on a set of measure zero for yij = ^'.J^ ^ , e.g. reference vectors which 

are "relative positions", with freely chosen "positions" yi G R 3 . In Section 5.2 we show that 
state synchronization with (a) is ill-conditioned. Simulations tend to indicate that this is 
characteristic of reference vectors defined as "relative positions". 

(ii) If k = 4 agents are not connected as in Theorem 4(c), then generically there is a con- 
tinuum of relative states corresponding to C Q \C S . Indeed: either one agent has a single link 
and therefore the necessary condition of Theorem 1 is not satisfied, or agents are connected 
in square like on Fig. 2. In the latter case, output synchronization means R\ = QiQ 2 G 
stabd/12), R 2 = Q 2 Ql G stab(y 23 ), R 3 = QiQj G stab(yi 4 ), R 4 = QaQI G stab(2/ 34 ) with 
consistency relation R\R 2 = R 3 R^. For generic y^, the first set of constraints allows R±R 2 
and -R3-R4 to span two 2- dimensional manifolds in ?>- dimensional 5*0(3) ; their intersection 
is generically 1- dimensional. Thus relative states (R±,R 2 ,R 3 ) characterizing C a span a 
one-dimensional manifold (whereas C s corresponds to R± = R 2 = R 3 = I 3 ). 

Example 5 The general theorems of Section 4-1 relate to Theorem 4 as follows. 

• For the setting of Theorem 4( c ); rankL 9 < j^E = 6 < nik — 1) = 9, so the condition of 
Theorem 2(a) is not satisfied. One checks that the condition of Theorem 2(b) holds. 

• For the setting of Theorem 4(b), condition of Theorem 2(a) is a fortiori still not satisfied, 
but the condition of Theorem 2(b) holds and ensures that C s is locally asymptotically stable, 
despite C Q 7^ C s . In fact C Q has two separated sets, one of which is C s . 

• Take, as an example setting of Theorem 4( a ), the situation of Example 1 but with refer- 
ence vectors y 12 = y\ 3 = (1,0, 0) T and y 23 = (0,1, 0) T . From Theorem 4( a ) an d Proposi- 
tion 2, C s = C is locally asymptotically stable. However, the condition of Theorem 2(b) is 



18 



Figure 3: Non-rigid graphs for which (A.) C a = C s from Corollary 1 and (B.) C Q ^ C s . 

NOT satisfied, so convergence cannot be exponential. 

The condition of Theorem 2(a) appears to be too strong for small networks. 

Theorem 4(c) gives a generic condition for C Q = C s on a small graph. The following 
algorithm exploits this condition by collapsing repeatedly subnetworks having the graph 
structure as in 4(c) to single nodes. If the resulting graphs consists of a single point, then 
generic reference vectors yield state synchronization. 

Algorithm 

• As long as there are i,l,m,pEV with (i,l),(i,m),(p,l),{p,m),(m,l) G £ 

(a) Collapse i,l,m,p to a single vertex % -- the resulting graph can 
have several edges connecting two vertices; 

(b) As long as there are vertices i,j connected by more than one edge: 

Collapse i,j into a single vertex i; 

• Output the resulting graph; 

Corollary 1 If the graph analysis algorithm above reduces Q(V,£) to a single vertex, then 
C D = C s for a generic set of reference vectors y^. 

Proof: If agents i, s are constrained to Qi = Q s for a network state in C , then all their 
links with the rest of the network are constraints on the same Q G 50(3). Thus we can 
collapse i, s into one agent and attribute it all edges incident from i or s. In addition, if 
agent j was connected to both i and s, then there is now a double-constraint Ryij = yij 
and Ryi S = yi s on R = QjQj = Q s Qf G 50(3); this generically requires R = J 3 , so j also 
belongs to the agents with orientations identical to Qi at C a . Our algorithm consecutively 
collapses agents in this way, according to (i) Theorem 4(c) on 4-agent groups, assuming 
generic y^, and (ii) the double-constraint implication just recalled. If this procedure can 
be applied until one agent remains, then C Q = C s . □ 

Remark 6 The condition expressed by our graph analysis algorithm differs from planar 
formation rigidity, see e.g. [11 J for a recent study. With 3 agents, rigidity is a weaker 
condition as it holds for graphs satisfying the condition of Theorem 4(b). With 4 agents, a 
graph satisfying Theorem 4(c) is minimal both for (local) rigidity and for C a = C s . With 
more agents, our objective becomes easier to achieve: e.g. the graph on Fig. 3. A ensures 
C Q = C s but is clearly not rigid. 
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5.2 Error sensitivity 

Achieving C s = C Q with less agents for coplanar y^ (Theorem 4(a)) than for linearly inde- 
pendent yij (Theorem 4(b)) may seem surprising. In fact, the two solutions (Q2Q1, Q3Q1) G 
{(/3,/3), (i?2*,-R3*)} for output synchronization in the proof of Theorem 4(b), correspond 
to the intersection of circles {R2II23 '■ R2 G stab(yi 2 )} and {R 3 y23 : R3 G stab(yi 3 )} on 
5 2 . In the coplanar case, the circles are tangent: this leaves (I 3 ,I 3 ) as only solution, but 
also indicates a possible loss of hyperbolicity of the dynamical system, which can have 
undesirable consequences. The following formalizes this on basis of error sensitivity. In 
Section 5.3, simulations further illustrate the effect on the convergence speed. 

We consider measurement errors replacing Qjyij by 

QjEijVij (18) 

with fixed perturbations E^ G {R G 50(3) : \\R — I 3 \\f < £ ^ 1}- The actual closed-loop 
evolution with these perturbed measurements, 

f t Qi = Qi Yl skew(Qj EijMijEjiQj) , % = 1, 2, . . . , k , (19) 

is the gradient of f _ e = J2(ij)es( l ~~ tr (Ql E ij M ij E JiQj))- A priori, Ey ^ Eji breaks 
the output map symmetry and Qi = Qj does not yield synchronization of the perturbed 
outputs. For clarity, we denote _ e (£'i2, . . . , Ej. (k-i)) the output consensus set with per- 
turbed outputs (18). For a large number of edges, / D _ e is strictly positive over (SO(n)) h 
for generic reference vectors and Ey, meaning that _ e = 0. We define 

k 

dist 2 (C A ,C B ) := min V] \\Q iA - Q iB \\ 2 F ■ 

{QlA,Q2A,—,QkA)&CA, 
{QlB,Q2B,---,QkB)&CB ~ 

Theorem 5 Consider k = 3 agents applying (19) on 50(3). Denote by E £ the set 
{(E 12 , . . . , E 32 ) : ||/ 3 - E^Wf < e Vz,j} C (50(3)) fc ( fc ~ 1 ). 

(a) In setting of Theorem 4( a )> there is a constant c > such that for any sufficiently 
small s > there is a non-zero measure set S £ C E e where C - e (E\2-, • • • , £32) is nonempty 
and dist(O s , C - e (E 12 , E 32 )) > c^fe for all (E 12 , • • • , E 32 ) G S e . 

Furthermore, for any e > there is also a non-zero measure set N e C E £ such that 
_ e (£ 12 , . . . , E 32 ) = for all (E 12 , . . . , E 32 ) G N £ . 

(b) In setting of Theorem 4(b), there is a constant c > such that, for any sufficiently 
small e > 0, C _ e (Ei 2 , . . . , E 32 ) is non-empty and dist(O s , _ e (£'i2, . . . , E 32 )) < ce for 
all (E 12 , . . . , E 32 ) G E e . 

Proof: The simple but tedious proof, constructing specific sets, is in the appendix. □ 

So in the setting of Theorem 4(a), the ratio between state synchronization error 
dist(O s , _ e ) and measurement error e can grow infinite for e — > 0, while in the set- 
ting of Theorem 4(b) this ratio remains finite; in this sense only the second situation is 
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robust. From geometric considerations (i.e. set intersections, see proof of Theorem 5), 
we also expect good robustness in the setting of Theorem 4(c). Simulations suggest that 
this distinction — y^ unrestricted robust, y^ restricted to relative positions not robust - 
carries over to Corollary 1. 

5.3 Numerical Simulations 

In this section we present some numerical simulations for 50(3). 

Figure 4 shows local evolution of f a and f s for 6 all-to-all connected agents implement- 
ing (6). The thick lines are for generic yij (denoted B, analogous to case(b) of Theorems 4 
and 5), thin lines are for y^ derived from generic positions yi G M 3 , that is y^ = ±^y~J*y.\\ 2 
(denoted A, analogous to case(a) of the theorems). Top plots are without measurement 
errors: for B, f a and f s exponentially decrease to machine precision, for A the decay is 
clearly non-exponential. Note that the conditions of Theorem 2 both fail for A, while 
they hold for B. Bottom plots repeat the simulation with random measurement errors 
\\Eij — Is\\f < £ = 0.01. As C _ e cannot be reached exactly, f _ e stabilizes between 10~ 4 
and 10 -3 for both cases. Case A now actually saturates at a significantly higher value of f s 
than B. This is reminiscent of Theorem 5: for case(b) the synchronization error is of order 
e, for case(a) it can be of order ^fe > e. The difference between subexponential conver- 
gence and saturation for A, becomes apparent on exceedingly long simulations (displayed 
on the small plots). The setting with y^ given by relative positions seems ill-behaved in 
practice. 

Figure 5 illustrates behavior with time- varying y^. Again we consider 6 all-to-all cou- 
pled agents (no measurement errors), with reference vectors derived from generic positions 
in 1R 3 , that is y^ = ±\\y.J y 3 .\\ 2 ■ The yi(t) vary quasi-periodically in time with frequencies in 
[0.05, ±] Hz. For e = 0.1 in (15), that is the thick lines, roughly exponential convergence 
to f — f s — is observed, with oscillations from time-variation. Note that this e is not 
exceedingly small w.r.t. frequencies of yiif). The value e — 10 (thin line) is too large for 
exploiting the time-varying setting: corresponding curves look more like the ones for fixed 
yij (thin lines on Figure 4), converging to very poorly, if at all. 

Finally, we illustrate the influence of graph Q on convergence. Figure 6 shows local 
evolution of f a and f s for 7 agents interconnected according to graphs A and B of Fig. 3, 
page 19. Both runs use the same generic fixed y^ (except for the missing link in graph B) 
and same initial conditions in the neighborhood of C s . The systems similarly converge to 
f — for both graphs. For graph A (thick lines), in accordance with Corollary 1, f s also 
converges to 0. For graph B (thin lines), f s converges to an arbitrary value, as C Q \ C s 
admits a continuum of relative states (Q2QJ > • • • , Q7Q1) 7^ (-^3, • • • , -^3)- 

6 Conclusion 

The present paper studies autonomous synchronization of orientation matrices Qi G SO(n) 
with coupling through outputs ^ states: information exchanged on each link is the expres- 
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Figure 4: Evolution of / Q _ e and f s for 6 all-to-all connected agents, random initial condi- 
tions: thick lines generic y^ (like case(b) of Th.4), thin lines y^ = ±^y~J*y.\\ 2 with generic 
yi (like case(a) of Th.4). Top plots without measurement errors, bottom plots with mea- 
surement errors. Far right plots: pursuing the simulation of case (a) to larger times. 
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Figure 5: Evolution of f a and f s for 6 all-to-all connected agents, random initial conditions 
and yij = ~±^Ezyh^~ 2 with time- varying yi{t). No measurement errors. Thick lines e = 0.1 in 
(15); thin lines e = 10. 

sion of a vector in respective body frames. A natural gradient law generalizing full-state 
coupling case [31] is proposed. Its convergence properties feature interesting specificities, 
even locally. These are studied in detail, with both necessary conditions and sufficient ones 
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Output disagreement — no meas.error State disagreement — no meas.error 




Figure 6: Evolution of f Q and f s for 7 agents with generic y^ G S 2 and initial conditions in 
the neighborhood of C s , with interconnection graph of Fig. 3. A (thick lines) and of Fig. 3.B 
(thin lines). No measurement errors. 

for state synchronization. An algorithm is proposed to find/build networks in which output 
synchronization implies state synchronization on SO (3) for generic reference vectors. The 
setting with reference vectors derived as relative positions in R 3 is shown to be non- robust. 
An analog system with states x,- t G R n is studied for comparison. 

We believe that this study can be both of practical interest — to design synchronizing 
networks with less than full-state exchange — and of theoretical interest — illustrating 
non-robust cases, highlighting further differences between R n and SO(n), and maybe mo- 
tivating a more detailed study of orbit intersections on homogeneous spaces and Lie groups. 
Regarding distributed systems, the non-trivial issues raised in this work motivate further 
study of observability-related theories for agents interacting in a network. 

7 Appendix 

Lemma 1 Consider an equilibrium Q := (Qi, ■ ■ ■ ,Qk) of (6). Define C r s := {(A, . . . , A) G 
(R nxn ) fc : A G R nxn } and F tj = QfM tj Qj. Build the symmetric matrix L" = (If/) G 

R knxkn L 99 = _ a .. F .. j Qr • ^ ■ ^ j99 = ^ Q . . p.. j,^ g = {(g^^ _ _ _ ^ Q fe Q*) . 

Q* G SO(n)} is a set of equilibria, and S is locally exponentially stable if and only if 
(I n <g) L") is positive definite on the subspace of vec(so(n) fc ) orthogonal to vec(CJ). 
(N.B.: the condition is necessary only for exponential stability.) 

Proof: At a critical point x the Hessian quadratic form Hf(x) : {T x Ai x T X M) — > R of a 
smooth function /: M. — > R on a smooth manifold A4 is fully defined by Hf(x)(v,v) := 
4j?(f o 7)(0) for any smooth curve 7(t) on M. with 7(0) = x and ^7|o — v , see [10]. On 
a Riemannian manifold one can define at critical points 4 the self-adjoint Hessian operator 
Hy(x): T x Ai — > T X M. by (v,Ej(x)v) = Hf(x)(v,v) for v G T^A^ with (•, •) the Riemannian 
metric. 

4 For other points one would have to use the Riemannian connection. 
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Take 7: R -> (SO(n)) fc with 7 (0) = and | 7 (0) = (QA,...,^) =: 

Qfi G TQ(SO(n)) k , with 0* G so(n); also denote Q = (fij, . . . , Q k ) G R nxnk . From (5), 

H fo (Q){QQ,QQ) = ti{pLL 9 m T ) = (vec(fi T )) T (J n ® L")(vec(ft T )) . 

Right-invariance of (6) implies that S* is a set of equilibria on which / G = f (Q) is constant. 
Moreover, curves with tangent vector in Q(so(n) k nCl) = {(Qiflb, ■ ■ ■ , Qk^b) '■ G so(n)} 
at any point Q = (Qi, . . . , Qk) G (SO(n)) k follow invariance directions of the dynamics. In 
particular, S is a submanifold generated by all of these curves that pass through Q. The 
condition on (I n <g> L") ensures that H fo (Q)(QQ, QQ) = implies Qtt G Q(so(n) k n C s r ) = 

As in the proof of Theorem 3 we project the dynamics onto the reductive homogeneous 
space Ai = (SO(n)) k /C s and denote by ir: (SO(n)) k —> M. the canonical projection. The 
function f Q is constant on the fibers of ir and induces a smooth function f on Ai which 
satisfies f = fo ^- Since (6) is invariant under right multiplication of all Qi by the same 
Q G SO(n), it induces a vector field X on Ai. The set S is collapsed by ir onto a single 
equilibrium s6M and it is locally exponentially stable under (6) iff s is locally exponen- 
tially stable under X. We equip Ai with the induced normal Riemannian metric (•, -)m and 
denote by f)(Q), Q G (SO(n)) k , the canonical horizontal distribution on (SO(n)) k , i.e. the 
orthogonal complement of the vertical distribution X>(Q) := ker(Tg7r) = Q(so(n) k D CJ) 
w.r.t. (•, -)m- By the invariance of f under the action of C s on (SO(n)) k by right multipli- 
cation, the vertical part (grad/ (Q))° G tj(Q) of grad/ (Q) is zero. Recall that for the nor- 
mal metric, (V, W)( S o(n)) k = (T$k(V), Tq1t(W))m for all Q G (SO(n)) k , V, W G fj(Q) [8]. 
Some calculations, using properties of the normal metric and the definition of the gradient, 
show that the projected vector field X(tt(Q)) is actually grad^ fo{^{Q)), the gradient of 
f w.r.t. the normal metric. So grad/ Q projects to a gradient system on Ai. Therefore, 
the equilibrium point s G Ai is locally exponentially stable if and only if the Hessian form 
of f at s is strictly positive definite (this follows directly from standard results on expo- 
nential stability and the fact that the Hessian operator at a critical point corresponds 
to linearization of the vector field in exponential charts). Simply using the definition of 
the Hessian form and lifting curves from Ai to SO(n) k horizontally, we see that 
is strictly positive definite at s if and only if Hf o is strictly positive definite on t)(Q). But 
from the above, Hf o is strictly positive definite on t)(Q) iff (J n (8> L") is positive definite 
on the subspace of vec(so(n) k ) orthogonal to vec(CJ). □ 

The following result is used to prove Theorem 5. Let h: SO(3) — > S 2 given by h(Q) = 
Q T y for y G S 2 and denote h" 1 ^) = {Q : Q T y = p} for p G S 2 . 



Proposition 5 If y, z G S are linearly independent, then there is a neighborhood U 0/J3 
in SO (3) such that for all p G S 2 , R G U, S G 5*0(3), the set h~ l (p) r\Rstab(z)S contains 
at most one element. 

Proof: Note that Rst&b(z)R T = stab(ite). One checks that h~ 1 (p)S T R T = h~\RSp). 
Thus (/i _1 (p) fl Rstab(z)S) = (h~ 1 (RSp) D stab(-Rz)) RS. Given linearly independent z, y, 
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there is a neighborhood U of J3 in SO (3) such that y and v := Rz remain linearly inde- 
pendent for all R GW. This linear independence implies that equation X T y = RSp has at 
most one solution X £ stab(i>) = stab(Rz) for all p £ S 2 . □ 

Proof of Theorem 5: C _ e is invariant under right multiplication of the states by a 
common constant, so (to simplify notations) we can characterize C _ e by first assuming 
Qi = I3 and at the end reintroducing a common right-multiplication. A (hypothetical) 
state (I 3 , Q2, Q 3 ) in the perturbed output consensus set C „ e corresponds to Q^EyiDvi — 
#212/12, Q3-E132/13 = #3i2/i3 and Q 3 Q 2 r #232/23 = #322/23- The first two constraints yield 
Q2 £ #21 stab (2/1 2) #12 and Q3 £ #31 stab (2/13 )-E'J 3 . The output of agent 2 for j/23 is then 
contained in the circle 

Z' 2 = {E 12 QE^E 23 y 23 : Q £ stab(j/ 12 )} = {y £ S 2 : 2/f 2 # 1 T 2 2/ = 2/^1 #232/23}. 

Similarly, output of agent 3 for j/23 is contained in the circle 

Z' 3 = {E 13 QElE 32 y 23 : Q £ stab(j/ 13 )} = {i/eS 2 : y^E^y = y^E^E^} . 

For a state in C _ e , these two outputs must be equal, so they must belong to Z' 2 H Z 3 . 
Since the pairs 2/12,2/23 and 2/13,2/23 are each linearly independent, each point p £ Z' 2 D Z 3 
corresponds to exactly one state of the form (I3, Q 2 , Q3) £ C _ e for sufficiently small 
perturbations \\Eij — Is\\f, by Proposition 5 (existence of the state is ensured by definition 
of Z' 2 and Z' 3 ). Moreover, 2/12,2/13 linearly independent ensures that Z' 2 D can contain 
at most two points for sufficiently small \\Eij — I 3 \\f- Let us introduce G 3 := Ej 1 Q 3 Ei 3 , 
G2 '■= Ef 3 E 12 E 2l Q 2 Ei 3 , y 12 = Ef 3 Ei 2 y 12 , y 23 = E 31 E 32 y 23 and E A = Ef 3 E 12 E 21 E 23 E 32 E 31 . 
A state being in C _ e then requires (G2,G 3 ) £ stab(j7i2) x stab(2/i3). Further, define 

Z 2 := E± 3 Z 2 = {y £ S 2 : j/^y = 2^ (#^23)} , Z 3 := E\ 3 Z' 3 = {y £ S 2 : 2/f 3 2/ = J/fsfe} • 

Again, there is a bijection between points in Z2 D Z3 and states in C _ e with Qi = / 3 - We 
also define the balls 

B 2 := {y £ S 2 : yj 2 y > yJ 2 (E A y 23 )} , B 3 := {y £ S 2 : y[ 3 y > y^s} 

whose boundaries are Z 2 resp. Z 3 . Furthermore, we denote by Z 2o , Z 3o the circles and 
B20, B 3o the balls for the unperturbed case, i.e. for = I 3 . We restrict ourselves to 
the case < 2/122/13 < 2/122/23 < 2/132/23; the proof for other cases is strictly analogous. In 
addition, we redefine d := J 12 ? Vl \ and y 23 = Jr^4^< ■ 

From invariance of the Frobenius distance, ||/ 3 — J2II.F < \\h — Ji\\f + \\h — J\ J%\\F for 
any Ji, J2 £ iSO^). This allows to show that for (#12, • • • , #32) £ #e, 

HQ — Je\\f < m£: if Je is the product of Q with m error rotations Eij (20) 

arbitrarily split into mi left and m2 right multiplications, mi + m2 = m. Also, xf (J 3 — 
Ji)%2 < ||-^3 — ^i||f f° r J\ £ SO(n), Xi,x 2 £ S* 2 . 
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(a) For Theorem 4(a), circles Z 2o and Z 3o intersect at a single tangency point 1/23; then 
(in our restricted case, else balls must be redefined) the same is true for balls B 2o and B 3o . 
First consider the second assertion of the theorem. For any e > we can find a /3 > 
and an open subset N £ C E £ of (5"0(3)) 6 such that 2/^2/23 > 2/132/23, 2/12-^2/23 > 2/122/23 + 
and 2/122/12 > 1 — f for all perturbations (£'12, • • • , -E32) G AT e . Then Z 2 and Z 3 are in the 
interior of B 2o and i?3 respectively, so they cannot intersect and C _ e = 0. 

Now consider the first assertion. For any e > sufficiently small, we can find an a > 
and an open subset N £ C i? e of (S*0(3)) 6 such that 

(*) 2/132/23 < 2/13^2/23 - oe, 

(ii) ?/r22/23 > yJ 2 E A y 23 + and 

("«) II2/23 - I/23 1 1 2 < fe 

for perturbations (£7i2, • • • , -B32) £ ^e- 

From (i), E A y 23 , on the boundary of B2, also belongs to interior of £> 3 , so i?2 H B 3 has 
non-zero measure and Z 2 D Z3 contains two points, corresponding to two separated sets in 
C _ e . Take a state (Q*, Q 2 , Q* 3 ) G C _ e with corresponding (G 2 *,G 3 *) and = G 3 *y 23 = 
G 2 *E A y 23 . Then (ii) implies j7f 2 j/* = yl 2 G 3 *y 23 = y[ 2 E A y 23 < yj 2 y 23 - as so j7f 2 (/ 3 - 
G 3 *)y 2 3 > and using (iii) we get y( 2 {I 3 - G 3 *)y 23 > fe. Writing G 3 „ and vectors 
in orthonormal basis (ei,j/i3,ei x z/ 13 ) shows that this requires G3* to be a rotation of 
angle 0(e) with (1 — cos(8(e))) > C\E for some constant C\ independent of e; then there is a 
constant C2 > independent of e such that 0(e) > c 2 y/e. Thus we can find, for e sufficiently 
small, a constant C3 > independent of e such that \\I 3 — G 3 *\\p = \/8sin(^-) > c 3 y/e. 
Using property (20) and 

dist(C„ C _ e ) = min (||g - h\\ 2 F + \\Q - Q 2 f F + \\Q - Q 3 \\ 2 F ) 1/2 

Q6SO(n) 

(l3,Q2,Q 3 )eC - e 

>j= min \\Q - I 3 \\ F + \\Q - Q 2 \\ F + \\Q - Q 3 \\ F 

(/3,Q2,Q3)eC* _ e 

min ||Q-/ 3 |^+ WQ-QsWf > \ min ||J 3 - Q 3 \\ F 

(^3,Q2,Q3)6C -e 

we get, for sufficiently small e, a constant C4 > independent of e such that 
dist(C s , C . e ) > - Q 3 4 F > IU3-G3«lbHIG 3 .-Q3.1| f > _y /3 _ Gs+ || F _ 2£ > ^ 

(b) Denote by 9 the angle between j/ 2 3 and y\ 2 i.e. 2/122/23 = cos(0), with G (0, |] since 
we consider 2/122/23 > 0. Define B„ = {1/ 6 : — 2/23 1 1 2 < «}• Theorem 4(b) requires Z 2o 
and Z 3o to intersect at two points. This requires their intersections to be transversal. As 
a consequence, there exist /?0)7o > such that a < 70 implies: (i) following the boundary 
of B a , one passes through 4 points p 2a , P3a, P2b, P3b in that order, with {p 2a ,P2b} C Z 2o and 
{P3a,P3b} C Z 3o ; and (ii) on each boundary arc [p»,Pj], there is at least one point 
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whose distances to Z 2o and Z 3o are both strictly larger than foot. Next, note that there can 
be no point belonging to Z 2 (resp. Z 3 ) whose distance to Z 2o (resp. to Z 3o ) is larger than 
sin(e/2) • Indeed, denote x G S 2 a point at distance 5 from Z 2o in K 3 , and denote <p G [0, 7r] 
the angle between vectors x and y 12 ; then | cos(6 l ) — cos(0)| > <5sin(^rp) > 5sin(|). At the 
same time, x E Z 2 requires y\ 2 [x — E A y 23 ) = hence 

I cos (0) - cos (</>)| = \yj 2 (x - y 23 ) - yj 2 (x - E A y 23 )\ 

= \yl 2 {h - E? 2 E 13 )x + y{ 2 {I 3 ~ E^E 23 )y 23 \ 

< llyialbdlis-^iallFlNb+ll/a-^i^allfll^lb) <4e. 

These two inequalities yield 5 < As/ sin(|). A similar argument applies for Z 3 . 

As Z 2 and Z 3 are thus confined to neighborhoods of Z 2o and Z 3o , take a < 70 and 
£ _ /3oasm(e/2) ^ Then the intersection points of Z 2 and Z3 with the boundary of B a 
must remain confined in their respective arcs delimited by the r\p isPj ], e.g. Z 2 has one 
intersection with the boundary of B a on the arc between r\p 2a>PZa \ and f[p 3i) ,p 2a ], similarly for 
the others; thus the order of intersections along the boundary — alternating between Z 2 
and Z 3 — is conserved. This implies that Z 2 and Z 3 have an intersection point y* G B Q , 
characterizing a state (7 3 , Q 2 *, Q3*) G C _ e associated to (G 2 *, G 3 *) G stab(yi2) x stab(yi 3 ) 
and = G 3 *y 23 = G 2 *E A y 23 . Then using property (20), 

dist(C s , C _ e ) = min (||g - J 3 ||| + \\Q - Q 2 f F + ||Q - Q 3 \\ 2 F ) 1/2 

Q&SO(n) 
(h,Q2,Qi)eC - e 

< min \\Q - I 3 \\f + \\Q - Qz\\f + \\Q ~ Qs\\f 

Qe50(n) 

(-f3,Q2,Q3)6Co-e 

< 11^3 ~~ Q2* || F + ~~ Q3* II F < ~~ ^2* \\f + ||-^3 ~~ G*3* \\f + 6£ . 

By the choice of a and B a 9 we have \\y* — y 23 \\ 2 < ^ ^nn^ ■ Moreover, as G 3 * is a 
rotation around y 13 , we have ||J 3 -(?a l ,||j' = ^\\y* -fjwh with r 3* = II2/13 x 2/23 lb ; similarly 
11-^3 - C2*||f = ^\\y* - E A y 23 \\ 2 with r 2* = 1 1 2/12 x (E A y 23 )\\ 2 . There exist ft, 71 > such 
that £ < 71 implies min(||f7i 2 x (E A y 23 )\\ 2 , \\y 13 x y 23 \\ 2 ) > \/2ft- Then 



,r r« \ Wy* — 2/23 1 1 2 + \\y* -E A y 23 \\ 2 2\\y* - 7/23II2 + 6£ 

dist(C s , C _ e ) < h 6£ < h 6£ < a 2 s 

Pi Pi 

for small enough s, with a 2 = 6 + (6 + ^ S in(g/2) ) nn it e - ^ 
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